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Abstract 

First principles microphysics models are essential to the design and analysis of high energy density physics experiments. 
Using experimental data to investigate the miderlying physics is also essential, particularly when simulations and 
experiments are not consistent with each other. This is a difficult task, due to the large number of physical models 
that play a role, and due to the complex (and as a result, noisy) nature of the experiments. This results in a large 
number of parameters that make any inference a daunting task; it is also very important to consistently treat both 
experimental and prior understanding of the problem. In this paper we present a Bayesian method that includes 
both these effects, and allows the inference of a set of modifiers which have been constructed to give information 
about microphysics models from experimental data. We pay particular attention to radiation transport models. The 
inference takes into account a large set of experimental parameters and an estimate of the prior knowledge through a 
modified function, which is minimised using an efficient genetic algorithm. Both factors play an essential role in our 
analysis. We find that although there is evidence of inaccuracies in off-line calculations of X ray drive intensity and 
Ge L shell absorption, modifications to radiation transport are unable to reconcile differences between ID HYDRA 
simulations and the experiment. 

Keywords: inertial confinement fusion, radiation hydrodynamic simulation, Bayesian inference, plasma opacity, 
uncertainty analysis, convergent ablator, national ignition facility, radiation transport 



1. Introduction 

In recent inertial confinement fusion (ICF) experiments 
performed at the national ignition facility (NIP) [T], sig- 
nificant differences between radiation-hydrodynamic sim- 
ulations and experimental data have been observed [2]. 
It is not clear whether these simulations are inaccurate, 
or that they neglect some important physical effect. It is 
quite challenging to investigate which aspects of physics 
models are causing discrepancies and should be improved, 
largely due to the complex, nonlinear dependance of ICF 
capsule evolution on a large number of underlying models. 

The complex nature of the experimental designs is an 
important source of the difficulties. There are a large 
number of experimental parameters that are only known 
with limited accuracy; variations in these parameters rep- 
resent a noise source in the experimental data that can 
reduce the significance of the experimental result. Since 
the physical models we aim to investigate are fairly well 
constrained by a large amount of previous work, it is im- 
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20 portant to take into account the relative significance of 

21 the experiment and of any previous work. The interplay 

22 between experimental and previous information is an es- 

23 sential ingredient in a reliable analysis, and is often ne- 

24 glected. Its inclusion requires a consistent treatment of 

25 all physical and experimental parameters; together there 

26 are far too many of these to treat directly, however they 

27 are too important to neglect completely. 

28 In this paper we present an analysis of experimental 

29 data taken from a single NIF shot, N110625. The aim 

30 is develop a method of investigating microphysics models 

31 taking into account many of the noise sources in the ex- 

32 periment, and prior work. We use an inference model that 

33 has been developed specifically to allow the large number 

34 of parameters to be dealt with in a consistent manner [3] . 

35 In this work we focus on inferring information about ra- 

36 diation transport in the ablator of an ICF capsule from 

37 time resolved data taken from radiography [?] . Radiation 

38 transport relies on several physics models which must be 

39 approximated to make a full capsule simulation tractable, 

40 and as a result are often considered to be potential sources 

41 of model inaccuracy. In this work existing microphysics 

42 tables are modified in physically motivated ways; these 

43 modifiers are interpreted as measures of the inaccuracies 

44 in the physics models, and their inferred values give in- 
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45 formation about the source of difficulties in describing 

46 experimental observations. 

47 2. Microphysics in Inertial Confinement Fusion 

48 ablators 

49 In a typical indirect drive ICF design a spherical plastic 

50 shell, filled with deuterium-tritium (DT) fuel, is bathed 

51 in X rays created by the interaction of laser light with 

52 a high Z hohlraum. The resulting ablation of the outer 

53 plastic produces a rocket action that implodes the shell, 

54 compressing the fuel until it undergoes thermonuclear fu- 

55 sion. The propagation and absorption of X ray energy 

56 in the plastic and fuel is an essential piece of describing 

57 the implosion that requires detailed models of microscopic 

58 physics. These physics issues are described by a suit of 

59 computer simulations which provide, for example, tables 

60 of radiative opacities which are taken as input by subse- 

61 quent radiation- hydrodynamic simulations [5 . 

62 Achieving ignition is a challenge and so the design of 

63 successful targets requires careful tuning of a large set of 

64 design parameters, based on the results of simulations [6 . 

65 This means that the fine details of microphysics models 

66 can be very significant; nevertheless, the important as- 

67 pects can be understood with relatively simple one dimen- 

68 sional models }7 . We will discuss the important aspects 

69 of microphysics models in these terms. 

70 At its peak, the X ray drive on the outer surface of the 

71 capsule has a brightness temperature of around 300eV. 

72 The majority of the energy of this field is in photon en- 

73 ergies that coincide with the K shell absorption edge in 

74 carbon (which accounts for ^-^50% by number of the plastic 

75 ablator) and so the model of this absorption feature plays 

76 an important role in determining energy deposition in the 

77 ablator. Higher photon energies are able to propagate all 

78 the way through the carbon, depositing their energy in 

79 the DT fuel. Heating of the fuel by these hard X rays 

80 has a detrimental effect on the implosion since, for the 

81 efficient adiabatic implosions driven by the NIF, the final 

82 density is in part determined by the initial temperature 

83 of the DT. Preheat by X rays reduces fuel compressibility 

84 and ultimately reduces the final convergence that can be 

85 acheived. An important player in this preheat is emis- 

86 sion from the M shell of the gold hohlraum wall, which 

87 produces an enhancement over the thermal specturm of 

88 photon energies > f .8KeV; in order to block these from 

89 reaching the fuel a dopant layer is buried in the ablator. 

90 In this work we consider germanium doped ablators, in 

91 which case absorption by the Ge L shell aligns with the 

92 Au M shell emission and prevents preheat of the fuel. 

93 In reality, the growth of 2 and 3 dimensional instabili- 

94 ties also plays a very important role in determining the 

95 implosion efficiency. In severe cases these can be much 

96 more important than the ID considerations that we have 

97 described. 



98 These two aspects of radiation transport, namely ab- 

99 sorption of the drive field and preheat of the fuel, clearly 

100 depend on models of the generation of the drive spectrum 

101 and of the absorption in carbon and germanium at a large 

102 range of conditions (10-200 eV, I-IO g/cc). They also have 

103 direct consequences for the dynamics of the implosion. A 

104 simple rocket model for the inwards acceleration of the 

105 ablator [7] shows that the velocity and remaining ablator 

106 mass are directly related to velocity of the ablated mate- 

107 rial, and therefore the absorption of drive radiation. The 

108 density of the fuel at a given time is related to the preheat. 

109 Measurements of these three quantities, as described in 
no section [4j can therefore provide information about under- 

111 lying radiation transport physics models. The complexity 

112 of ICF experiments and radiation-hydrodynamic simula- 

113 tions means that extracting this information is a chal- 

114 lenging data analysis problem; we describe a method of 

115 performing this analysis in the next section. 

116 3. Bayesian analysis of ICF experiments 

117 The relationship between physical models, which them- 

118 selves are very complex, and the data is approximated by 

119 radiation-hydrodynamic simulations which may not be 

120 well behaved enough to allow the use of computational 

121 inversion techniques f8^ or fitting techniques [9l [10] ; the 

122 large number of physical models that control the evolution 

123 of an ICF target also presents a problem for these meth- 

124 ods. The complex nature of the experiments also means 

125 that there are a large number of experimental parameters; 

126 although these are often constrained by target metrology 

127 and design tolerances, their large number makes them 

128 a significant source of noise in simulations [llj . Dealing 

129 with the very large space of physical and experimental pa- 
no rameters is an important challenge to a consistent analysis 

131 of ICF data. The usual methods of reducing the number 

132 of parameters, for example by Monte-Carlo sampling (see, 

133 for example, fl2!), are prohibitively expensive, and simply 

134 neglecting parameters will lead to misleading results. 

135 In l3] we have developed an inference method that al- 
oe lows these problems to be addressed. The approach is to 

137 separate out those parameters that are known to affect 

138 radiation-hydrodynamic simulations but are not of direct 

139 interest to the investigation of microphysics models; these 

140 are defined as 'nuisance parameters'. Typically these pa- 

141 rameters refer to experimental variables which have a 

142 known probability distribution, for example a target di- 

143 mension that has been measured with some error bar. 

144 The probability distributions of all nuisance parameters 

145 are mapped onto the output of radiation-hydrodynamic 

146 simulations; as a result the simulation output can be con- 

147 sidered as being probabilistic. In our model we assume a 

148 linear response to nuisance parameters, resulting in an an- 

149 alytic expression for the probability distribution of simu- 

150 lation outputs (the likelihood). Parameters that are phys- 

151 ically interesting (and therefore will be inferred from ex- 



2 



152 perimental data), such as microphysics models, are kept 

153 separate from the nuisance parameters aUowing their re- 

154 lationship with experimental data to be described using 

155 the full complexity of the simulation code. 

156 The inference model we have outlined is based on the 

157 maximum a posteriori (MAP) estimate; that is, the most 

158 probable values of all parameters of interest when the 

159 experimental data and prior have been taken into account. 

160 In our analysis these values are found by minimising the 

161 function [S] 



exp,t 



- id,,p - drniO)) 13' I3{d,,p - d^{e)) 

+ iln(|A,||a^a|)-lnP(^?) (1) 

with respect to the vector of interesting parameters 6. 
In the above expression, dm{0) is the vector of simula- 
tion outputs for given interesting parameters and nomi- 
nal values of the nuisance parameters, d^,j.p is the vector 
of experimental data, P{0) is the prior distribution of in- 
teresting parameters (discussed below) and the matrices 
a and /3 satisfy the equations 



A; 



162 These matrices summarise the effect of nuisance param- 

163 eters on our analysis; A^xp and A^ are the covariance 

164 matrices of the experimental measurement and nuisance 

165 parameters, respectively, and A is the linear response of 

166 the simulation to nuisance parameters 77: A^j = ^'^g^^'^^ . 

167 Equation ([T]) takes the form of a modified function. 

168 The first term on the right hand side is the usual anal- 
og ysis, and the second can be interpreted as a loss of infor- 

170 mation from the experiment due to nuisance parameters. 

171 The third is a normalisation factor. The final term ex- 

172 presses the contribution from prior work on the values of 

173 the interesting parameters. In our application we inter- 

174 pret this term as an estimated error bar on the physical 

175 models we aim to investigate, reflecting previous work to 

176 validate them. The inclusion of this prior information 

177 provides context for the experimental result, allowing in- 

178 ferences to be obtained from a single observation. In [Hj 

179 this was shown to play a very important role in the anal- 

180 ysis of NIF data. 

181 The summary of nuisance parameters in the matrix 

182 0^13 has reduced the number of variables we must con- 

183 sider to only the ones of direct interest. The resulting 

184 smoothing of the simulation output also means that the 

185 minimisation of equation ([T]) can be approached using 

186 standard numerical methods. In this work we use a ge- 

187 netic algorithm (GA) to efficiently perform the minimi- 

188 sation. The details of the genetic algorithm have been 



189 optimised to allow an efficient exploration of a large pa- 

190 rameter space; the sacrifice is that the algorithm is more 

191 likely to find local minima. In the case of the ICF data 

192 we will consider here, this is not expected to be an issue 

193 since the interplay between likelihhod and prior tends to 

194 produce a single minimum. In more complex cases this 

195 can be tested by using several random initialisations, or 

196 avoided by using a more robust algorithm. 



197 4. Application to NIF experimental data 

198 We aim to demonstrate the application of our Bayesian 

199 inference method to the investigation of microphysics 

200 models using NIF data. We use ID simulations of a cap- 

201 sule implosion performed using the HYDRA radiation- 

202 hydrodynamics code jl3j . Our investigation proceeds by 

203 defining a set of modifiers to the inputs of these simula- 

204 tions, and inferring the values of these modifiers. We are 

205 concentrating on physics issues in radiation transport and 

206 so our modifiers are to the X ray drive spectrum imping- 

207 ing on the capsule's outer surface (found from seperate 

208 models of the hohlraum) , and to relevant opacity models 

209 of the ablator material (taken from the TABOP opacity 

210 model). The use of modifiers, placed on the results of 

211 existing calculations, allows our inference results to be 

212 interpreted as implied inaccuracies in microphysics mod- 

213 els. We give details of our modifiers in table[lj In the case 

214 of the drive timing modifier, the prior error bar reflects 

215 the error bar on the DANTE instrument 14 , which gives 

216 a time-resolved measurement of the drive radiation tem- 

217 perature. This instrument has played an important role 
21s in the development of the separate hohlraum simulations 

219 which produce drive proflles for our capsule simulations. 

220 For all other modiflers, prior errors are estimated in order 

221 to reflect the expected accuracy of the underlying physi- 

222 cal models. All modifiers, with the exception of the drive 

223 timing, are dimensionless multipliers on existing models 

224 and so their 'nominal' (and therefore prior) values are 1; 

225 the drive timing has a nominal shift of ns. 

226 Experimental data are taken from a single NIF 'con- 

227 vergent ablator' shot, N110625. This experiment utilised 

228 a germanium doped capsule which was radiographed as 

229 it imploded giving a time- and space- resolved measure- 

230 nient of plasma density jH I15j . This then gives time- 

231 resolved data for the implosion velocity, mass of the ab- 

232 lator, and the pR product of the imploding fuel shell. 

233 We consider these three data points, taken at three times 

234 during rocket-like phase of the implosion, in our analysis. 

235 The use of implosion velocity and ablator mass, which 

236 diagnose the drive, along with the pR which is sensitive 

237 to preheat of the fuel, should allow the degeneracy of our 

238 modifier set (for example the drive intensity and C K 

239 shell) to be lifted. This is important since such degen- 

240 eracy results in a set of multiplier values that minimise 

241 equation ([T]); the inclusion of the pR data should select a 
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Modifier Name 


Description 


Expected Effect 


Prior Error 


Drive Intensity 


Multiplies intensity of 4*^ rise in X 
ray drive 


Increased drive results in increased ve- 
locity and decreased ablator remaining 
at given time 


±0.1 


Drive Timing 


Shifts the timing of the 4"^ rise of 
the X ray drive 


Earlier rise increases drive at given time 


±0.1 


Au M Shell 


Multiplies the intensity of the gold 
M shell component of X ray drive 
spectrum 


Increased M shell results in increased 
preheat and reduced pR at given time 


±0.2 


C K Edge 


Multiplies the opacity of the K shell 

absorption edge in carbon 


Increased absorption increases effective 
drive 


±0.1 


Ge L Edge 


Multiplies the opacity of the L shell 
absorption edge in germanium 


Increased absorption reduces preheat 


±0.1 



Table 1: Description of the modifiers placed on input physics models. The values of these modifiers are inferred from 
experimental data using the method described in the text, and are intended to give information about the accuracy 
of radiation transport models for NIF ablators. 



Modifier 


No Prior 


Including Prior 


Drive intensity 


0.57 


0.90 


Drive timing 


-0.45 ns 


0.01 ns 


Au M shell 


1.84 


0.97 


C K edge 


0.92 


1.0 


Gc L shell 


1.15 


1.16 



Table 2: Positions of the best fit to experimental data for 
NIF shot N110625. In both cases 29 nuisance parameters 
are included; comparison of the two sets of data measures 
the significance of the experimental data when compared 
to prior knowledge. 



242 single one of these values since it more fully refiects the 

243 physics of the problem. 

244 For the multipliers and experimental data described, 

245 our genetic algorithm is randomly initialised and procedes 

246 by automatically calling HYDRA. The nuisance param- 

247 eter modification 0^13 is calculated for the 29 physical 

248 dimensions, densities and material composition parame- 

249 ters of the capsule [5], which are assumed to be known 

250 with an error of 1%. We ran the GA for 25 generations 

251 with 92 members per generation, requiring up to 2300 

252 HYDRA simulations (the actual number is lower due to 

253 the optimisations made to the GA), equivalent to < 200 

254 CPU hours. In table [2] we give the position of the results 

255 for two cases; including and neglecting the prior, respec- 

256 tively. Since the position of the minimum of equation ([T]) 

257 is determined by the relative importance of the prior and 

258 experimental results, comparison of these two cases pro- 

259 vides information about the significance of the experiment 

260 in measuring radiation transport physics. 

261 In table [2] the fit given in the 'No Prior' column corre- 

262 sponds to a maximum likelihood (ML) analysis, in which 

263 the experimental data are the only source of information 

264 about the values of the modifiers. In this case, the results 



265 demonstrate that in order to fit the data all modifiers 

266 should be significantly different from their nominal values; 

267 this implies that microphysics models are in considerable 

268 error. Given the extensive work that has been under- 

269 taken on these models in the past, it is unlikely that this 

270 is truly the case. The previous work is taken into account 

271 in the 'Including Prior' column, and the large difference 

272 in results demonstrates the importance of including prior 

273 knowledge. In that case (corresponding to the MAP re- 

274 suit) all modifiers are much closer to their nominal values. 

275 The noise in the experiment makes the observed data in- 

276 sensitive to the details of radiation transport; only the 

277 overall drive and Ge absorption are significantly modified 

278 from their prior values. Our results suggest that the calcu- 

279 lated drive intensity is too high, consistent with previous 

280 work on IGF data, and that the calculated absorption by 

281 the germanium dopant layer is too low. 



282 Comparison of the best fits to experiment, found using 

283 the two inference methods (neglecting and including the 

284 prior) , allows us to measure the ability of inaccuracies in 

285 radiation transport to explain problems with modelling of 

286 the experiment. The quality of the inferred fits to experi- 

287 mentally inferred implosion velocity, ablator fraction, and 

288 line density are shown in figure [IJa),(b) and (c) respec- 

289 tively. In these figures, experimental data as a function 

290 of time are shown in blue, and simulation results using 

291 modifier values from table [2] are plotted in red (no prior) 

292 and black (prior included). The ML analysis, neglecting 

293 the prior, gives a reasonable qualitative fit to the data, 

294 but does not match within all error bars. The MAP re- 

295 suit is much closer to an unmodified simulation and gives 

296 a poorer agreement with experiment. The inability of ei- 

297 ther approach to give a good match to the data suggests 

298 that discrepancies between simulations and experiments 

299 are not solely due to issues with radiation transport. 
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300 5. Discussion and Conclusions 



a 

^ 240 




21.2 21.4 21.6 21.8 22.0 

Time (ns) 
(a) Implosion velocity 



0.22 
0.20 
0.18 
0.16 




21.2 21.4 21.6 21.8 22.0 

Time (ns) 

(b) Fraction of ablator remaining 



0.08 




21.2 21.4 21.6 21.8 22.0 

Time (ns) 
(c) Fuel pR 

Figure 1: Best fits to experimental data, corresponding 
to HYDRA simulations using modifiers given in table [2j 
Experimental data are shown in blue, and inference re- 
sults neglecting and including prior knowledge are shown 
in red and black, respectively. 



301 The work presented in this paper demonstrates a 

302 method for inferring information about first principles 

303 physics models from ICF data. The inference model we 

304 use allows the inclusion of a large number of nuisance 

305 parameters; these play an important role in determin- 

306 ing the information in the experimental result. This is 

307 essential when comparing experimental results with the 

308 results of previous work, which is often the case in high 

309 energy density physics. Although we focus here on ra- 

310 diation transport in ICF ablators, the issues we discuss 

311 are important in many of the experiments performed in 

312 high energy density physics, and the inference method we 

313 describe is easily applicable to any of these. 

314 The main result of this work is that prior knowledge 

315 about microphysics plays a very important role. Includ- 

316 ing this in a consistent manner allows meaningful informa- 

317 tion to be extracted from data, so that when data imply 

318 a modification to physics models the result truly reflects 

319 the state of the art. We have also shown that the complex 

320 nature of ICF experiments means that the neglect of nui- 

321 sance parameters and/or prior information in a simple 

322 or maximum likelihood analysis will give misleading re- 

323 suits. In this work 29 parameters have been varied by 1% 

324 in order to produce the information loss due to nuisance 

325 parameters; for the well characterised targets used at the 

326 NIF certain capsule dimensions are known to a much bet- 

327 ter level than this, however prior knowledge will play an 

328 essential role regardless. 

329 Once these factors are accounted for, there is evidence 

330 that both the overall X ray drive and the absorption of the 

331 germanium L shell are inaccurate. This could serve to fo- 

332 cus subsequent investigation of the underlying models (for 

333 example further inferences of inaccuracies in charge state 

334 balance) , however the poor agreement between the cur- 

335 rent best fit and the experimental data shows that issues 

336 with radiation transport cannot explain discrepancies be- 

337 twcen the details of ICF implosions and simulations. It is 

338 important to note that inferences based on an incomplete 

339 set of modifiers, which appears to be the case here, may 

340 never give a good fit to data. Until a good fit is found the 

341 physical meaning of multipliers is limited, and inferred 

342 values should be treated accordingly. 

343 The method used here has been specifically designed 

344 so that an analysis with a large enough set of modifiers is 

345 feasible. Cases with 1-2 orders of magnitude more eval- 

346 nations of ^re possible with a fairly modest computa- 

347 tional requirement, and the number of nuisance parame- 

348 ters can be increased in our linear model with almost no 

349 numerical overhead. 

350 Genetic algorithms have been previously used for 

351 HEDP applications, with good results [T5HIH]- In par- 

352 ticular, there is interest in using multi-objective genetic 

353 algorithms to consider several data sets simultaneously 

354 (typically 3 or 4). For the 9 data points we consider here. 
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355 and the even larger sets we aim to use, such muhiobjective 404 [7] S Atzeni and J Meyer-ter Vehn. The Physics of In- 

356 approaches would be difficult. Our single objective mod- 405 ertial Fusion, volume 125 of Internation Series of 

357 ified approach is in effect a linear scalarisation of the 406 Monographs on Physics. Oxford Science, 2004. 

358 multiobjective problem and allows much larger datasets 

359 to be considered. The trade off is that a single solution ™' t^l KM Hanson and GS Cunningham. Posterio sampling 

360 is found where multiobjective methods give several can- ^i^h unproved efficiency. Proceedings of the SPIE, 

361 didates; our careful treatment of the error bars on each ™ 3338:371-382, 1998. 

362 data point serves to justify our choice of scalarisation. [gj j g^^j^^^ ^^^^^^ Mitchell, and HP Wynn. 

363 It has been previously noted that the linear model we Besign and analysis of computer experiments. Sta- 

364 employ is not justified for ICF targets, since they have ^^^f-^^^i Science, 4(4):409-435, 1989. 

365 been highly tuned to operate at peak performance. The 

366 advantages of the analytic expression ([I]) are great, and 413 [10] Marc C. Kennedy and Anthony O'Hagan. Bayesian 

367 so the authors aim to develop an analytic model that is 414 calibration of computer models. Journal of the Royal 

368 more suited to ICF data. The linear model does, however, 415 Statistical Society: Series B (Statistical Methodol- 

369 capture the essence of the problem; that complex exper- 416 ogy), 63(3):425-464, 2001. 

370 iments produce less significant results when compared to 

371 existing knowledge. - [11] D. Clark. Capsule modehng of ConAbl N101220. 

™ • , r iu J 11 j-u • lis Presented at NIC Workshop, April 18 2011. 

372 Ihe Bayesian nature 01 our method allows the consis- ^ 

373 tent analysis of all available data, either by evolving the [^2] Byron P. Roe. Statistical errors in monte carlo csti- 

374 prior knowledge as more data becomes available or by j^ates of systematic errors. Nuclear Instruments and 

375 including all data in a single analysis; the different sets Methods in Physics Research Section A: Accelerators, 

376 of data do not need to be from the same experiment, or Spectrometers, Detectors and Associated Equipment, 

377 even ones of the same design. These extensions will form 570(1)T59 - 164 2007 

378 a important part of our further work. Finally, the compu- 

379 tational methods we have presented are suitable for both 424 [13] M. M. Marinak, G. D. Kerbel, N. A. Gentile, Et 

380 experimental design and discovery purposes, and we aim 425 al. Three-dimensional HYDRA simulations of na- 

381 to develop this application. 426 tional ignition facility targets. Physics of Plasmas, 

427 8(5):2275-2280, 2001. 
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